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ABSTRACT 


The  established  concept  of  temperature  determination  from  dual  spectral  band 
radiometric  measurement  is  analyzed  for  the  case  that  the  measurement  device  is  sensi* 
tive  to  total  photon  emittance  rather  than  radiant  emittance.  A  temperature  dependent 
ratio  of  measurements  is  calculated  for  two  distinct  spectral  bands  using  a  black  body 
distribution.  It  is  shown  that,  if  the  bands  are  not  self  contained,  then  there  exists  a 
one*to*one  correspondence  between  temperature  and  ratio. 

A  prototype  algorithm  is  proposed  and  tested  which  demonstrates  the  feasibility  of 
extracting  both  the  temperature  and  range  by  using  three  distinct  spectral  band  meas¬ 
urements.  The  model  assumes  nonhomogeneous,  wavelength-dependent  atmospheric 
attenuation.  The  target  and  sensor  systems  each  have  an  arbitrary  location  in  a  vertical 
plane  relative  to  the  earth.  In  the  computer  simulations,  LOWTRAN  data  was  used. 
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I.  INTRODUCTION 


A.  REVIEW  OF  PRIOR  WORK 

In  the  photon  detection  process  electromagnetic  radiation  causes  three  kinds  of 
mechanism  such  as  photon  efTects,  thermal  effects,  and  wave  interaction  effects 
[Ref.  1].  The  focus  of  the  remaining  discussion  will  be  confined  to  the  Arst  two  effects 
which  cover  the  majority  of  applications  involving  the  detection  of  thermal  radiation. 
One  common  photon  effect  is  photoconductivity.  Photoconductivity  is  observed  in 
semiconductors  and  characterized  by  photon  energy  and  energy  gap.  One  common  de¬ 
vice  which  uses  a  thermal  effect  is  the  bolometer.  Bolometers  can  be  made  of  any  ma¬ 
terial  which  has  a  temperature-dependent  resistance.  This  effect  is  characterized  by  the 
change  of  resistance  from  the  heating  effect  of  incident  radiation.  Therefore,  Plank's 
equation  for  radiant  emittance  can  only  be  applied  directly  to  thermal  effects.  Although 
some  devices  have  a  detection  mechanism  based  on  the  thermal  heating  effects  due  to 
radiation,  this  is  not  true  of  semiconductor-based  detectors.  Semiconductors  devices  are 
characterized  by  an  energy  gap.  Only  photons  with  energies  greater  than  this  gap  can 
produce  electron  hole  pairs  which  contribute  to  the  detected  signal.  The  excess  photon 
energy  is  dissipated  as  heat  [Ref.  2]  and  therefore  does  not  contribute  to  the  electrical 
signal  generated.  Specifically  the  integrated  spectral  emittance  [Ref.  3],  i.e.,  photon  flux, 
is  a  more  appropriate  signal  indicator  than  the  radiant  emittance  for  measurements 
based  on  semiconductor  devices. 

R.B.  Johnson  and  E.E.  Branstetter  presented  a  numerical  method  for  the  integration 
of  Plank's  equation  [Ref.  4].  Johnson  derived  an  approximated  equation  for  effective 
radiant  sterance  [Ref.  5].  Normalization  of  spectral  bandwidth,  for  convenience  in 
analysis,  was  provided  by  F.E.  Nicodemus  [Ref.  6].  Based  on  Plank's  equation  for  ra- 
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diant  emittance,  M.H.  Herman  presented  the  technique  for  temperature  determination 
using  the  ratio  of  output  from  two  bands  (Ref  7].  On  the  other  hand,  analysis  of  the 
dual  spectral  band  temperature  measurement  for  devices  based  on  the  photon  effects  has 
not  been  developed. 

One  standard  method  for  the  determination  of  the  gray-body  emissivity  and  the 
temperature  of  an  object  is  based  on  two  separate  narrowband  radiometric  pow'er 
measurements.  In  combination  with  Plank's  radiation  equation  for  a  gray-body,  two 
narrowband  measurements  are  sufficient  to  produce  a  unique  value  for  both  temperature 
and  emissivity.  For  practical  reasons,  such  as  improvement  of  signal-to-noise  ratio, 
most  instrumentation  take  measurements  over  bands  which  are  too  wide  to  be  consid¬ 
ered  narrowband.  In  this  more  general  situation  the  question  of  uniqueness  of  temper¬ 
ature  requires  more  careful  consideration.  This  issue  has  been  recently  analyzed  for  the 
distribution  based  on  a  Plank's  spectral  radiant  emittance  [Ref  8].  The  authors  devel¬ 
oped  an  analysis  based  on  the  evaluation  of  the  ratio  (^)  of  two  measurements  of  the 
integrated  power  distribution  over  separate  bands.  Three  cases,  nonovcrlapping  bands, 
partially  overlapping  bands  and  totally  overlapping  bands,  were  considered.  The  analyt¬ 
ical  investigation  demonstrated  that  the  last  case  docs  not  always  generate  a  one-to-one 
correspondence  between  the  ratio  ^  and  the  temperature. 

For  the  passive  ranging  technique,  P.J.  Ovrevo  and  R.C.  Wood  presented  the 
method  using  absorption  property  of  the  atmosphere  [Ref  9],  and  J.R.  Jenness,  Jr.,  and 
F.J.  Shimukonis  proposed  the  principle  of  passive  ranging  using  the  ratio  of  the  signals 
in  narrow  spectral  bands  received  [Ref  10]. 

For  optical  properties  of  the  atmosphere,  R.A.  McClatchey  and  others  described 
effects  of  index  of  refraction  in  the  atmosphere  [Ref  11:  pp.  41).  It  was  found  that 
transmission  factors  are  affected  by  no  more  than  1  -  2  %  if  refraction  is  included  in  the 
calculation. 
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B.  OVERVIEW  OF  THE  THESIS 


Chapter  I  is  the  introduction  and  covers  both  a  review  of  relevant  literature  and  a 
thesis  overview.  In  the  interest  of  keeping  this  thesis  fairly  self-contained,  Chapter  II 
provides  a  brief  technical  description  of  essential  background  preliminaries.  In  Chapter 
III,  effects  due  to  atmospheric  attenuation  and  detector  responsivity  are  ignored.  Using 
dual  spectral  measurements,  it  is  shown  that  it  is  possible  to  uniquely  determine  a  tar¬ 
get's  temperature  if  the  detector  bands  are  not  100  %  overlapping.  Specifically,  the  ratio 
of  the  dual  spectral  band  measurements  is  shown  to  exhibit  a  mathematically  predictable 
one-to-one  correspondence  with  the  temperature,  as  long  as  the  bands  are  not  totally 
contained  within  the  other. 

In  Chapter  IV,  a  numerical  program  for  calculating  a  spectral  band  measurement 
without  the  assumption  of  constant  atmospheric  attenuation  and  detector  responsivity 
is  discussed  and  tested  In  the  process  of  this  development,  a  novel  definition  for  a 
"responsivity",  applicable  to  devices  sensitive  to  photon  emittance,  has  been  introduced. 
If  the  effects  of  atmospheric  attenuation  on  the  dual  spectral  band  measurement  arc  not 
equivalent,  then  the  ratio  scheme  for  temperature  determination,  discussed  in  Chapter 
III,  is  not  applicable.  It  is  precisely  under  these  conditions  that  by  using  three  distinct 
bands  both  the  temperature  and  the  range  can  be  determined. 

In  Chapter  V,  a  prototype  algorithm  is  discussed,  which  demonstrates  the  feasibility 
of  obtaining  both  the  temperature  and  the  range  of  a  target  using  three  distinct  spectral 
band  measurements.  Complications  introduced  by  not  assuming  constant  altitude  of  the 
target-to-sensor  trajectory  arc  addressed.  LOWTRAN  data  was  also  used  here  to  char¬ 
acterize  atmospheric  attenuation.  Conclusions  for  the  thesis  are  presented  in  Chapter 
VI.  In  order  to  not  obscure  the  main  points  of  the  thesis,  the  extensive  mathematical 
derivations  have  been  relegated  to  appendices.  Both  Appendix  A  and  Appendix  B  fit  this 
description.  Appendix  C,  Appendix  D  and  Appendix  E  have  been  created  for  those 
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readers  who  are  interested  in  an  algorithmic  description  of  programs  used  in  the  com¬ 
puter  simulation. 
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11.  BACKGROUND 


In  this  chapter,  an  alternate  approach  based  on  photon  emittance  rather  than  radi¬ 
ant  eniittance  is  applied.  This  kind  of  treatment  is  appropriate  for  photovaltaic  or 
photoconductive  semiconductor  devices  but  not  radiometers.  First,  a  few  basic  as¬ 
sumptions  need  to  be  presented.  The  discussion  starts  from  the  ideal  case  with  the  fol¬ 
lowing  preliminary  assumptions  [Ref.  8:  pp.  1256]: 

•  opaque  object 

•  Lambertian  surface 

•  object  in  thermal  equilibrium 

•  spectrally  constant  emissivity  (gray  body) 

•  photon  emittance  derived  from  Plank's  radiation  law 

•  unity  spectral  transmittance  of  intervening  media  between  the  object  and  measuring 
instrument. 

Plank's  radiation  law  gives  the  spectral  radiant  emittance,  which  is  given  by 
[Ref.  3:  pp.  35): 


IV, ^,D 


2nc{}.)hc^ 

A^expihclUl')  -  1] 


[W •  cm~^  •  finr^], 


(2.1) 


where  e(A)  is  the  spectral  emissivity,  h  is  Plank's  constant,  X  is  the  wavelength,  c  is  the 
speed  of  light,  k  is  Boltzmann's  constant,  and  T  is  the  object's  temperature.  In  this  sec¬ 
tion,  since  the  emissivity  is  assumed  to  be  spectrally  constant,  e(X)  is  denoted  as  e. 

The  photon  emittance  is  derived  from  (2.1)  dividing  by  hclX,  which  is  the  energy 
associated  with  one  photon  [Ref.  3:  pp.  38].  Therefore  photon  emittance  is  the  flux  re¬ 
lated  to  the  number  of  photons.  The  resultant  photon  emittance  equation  is  obtained 
as 
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6a  ^ 


2ncc 


A^xpiftc/UT)  -  1] 


Iphotons  •  5“’  •  cm“^ ' 


(2.2) 


The  integral  of  (2.2)  over  some  band  limits  gives  the  photon  emittance  in  the  bands. 
The  photon  emittance  within  the  spectral  band  limits  [/li,  /Ij]  is  given  by 


2(^2,  A,.  7)-  (2.3) 

•’^1 

From  the  assumptions,  emissivity  is  assumed  to  be  spectrally  constant.  Therefore 
photon  emittance  for  a  gray  body  is  simply  expressed  as  photon  emittance  for  a  black 
body  multiplied  by  a  constant  emissivity  and  given  by 

6a  “  fiCsa/  (2.4) 


From  (2.3),  it  follows  from  (2.4)  that 


where  and  QbbP',  T)  is  the  spectral  photon  emittance  for  a  blackbody  [e(^)  —  IJ. 

It  has  been  shown  that  the  temperature  of  an  object  can  be  determined  by  taking  a 
ratio  of  measured  radiant  emittance  of  an  object  in  each  of  two  spectral  bands. 
[Ref  8:  pp.  1256]  In  Chapter  II  and  Chapter  III,  this  concept  is  applied  in  a  similar 
manner  except  the  assumption  is  made  that  a  ratio  is  formed  from  the  measured  photon 
emittance.  Let  the  two  spectral  bands  be  denoted  as  [.1],  A,]  and  [A4,  A,] .  Let  the  ratio 
for  photon  emittance  be  defined  as 
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(2.6) 


^’2*  ^3»  ^ 


Qi^.4,h.r)  • 


Since  from  the  assumption,  emissivity  is  unifonn  over  the  spectral  band,  einissivity  is 
cancelled  out,  and  the  ratio  becomes  the  ratio  of  the  photon  emittances  associated  with 
a  black  body,  as  shoum  by 


^4*  ^ 


l2aa(^2t  ^it  7) 
Qbb(^-4> 


(2.7) 


This  implies  that  di,  is  a  function  of  the  object's  temperature  T  only  and  the  ratio  can 
indicate  the  object  temperature.  When  the  dual  spectraUband  photon  emittance  method 
is  used,  it  is  critical  to  know  that  the  ratio  and  the  temperature  have  a  one-to-one 
correspondence.  For  this  purposes,  three  general  cases  of  spectral  bands  which  cover 

all  possibilities  are  considered.  The  possible  cases  ere  [Ref  8:  pp.  1256] 

1.  non-overlapping  bands  (Figure  la), 

2.  partially  overlapping  bands  and  common  band  limits  (Figure  lb,c),  and 

3.  totally  contained  band  (Figure  Id). 

These  three  cases  are  presented  graphically  in  Figure  1.  In  the  following  two  chapters 
the  validity  of  temperature  determination  based  on  dual  spectral  photon  emittance  for 
the  ideal  case  is  discussed. 
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a.  nonoverlapping  case 


1  2  a 

b.  partially  overlappii 

4 

ig  case 

13  2  4 

c.  one  common  band  limit  case 

1,3  2  4 


d.  totally  contained  case 


13  4  2 

Figure  1.  graphical  representation  of  possible  band  combination:  three  cases 
nonoverlapping,  partially  overlapping,  and  totally  contained. 


III.  DUAL  SPECTRAL  BAND  MEASUREMENTS  FOR  TEMPERATURE 

DETERMINATION 


In  this  section  the  ratio  of  photon  emittance  defined  in  (2.6)  is  considered  for  the 
various  cases  discussed  at  the  end  of  Chapter  II.  In  order  to  validate  the  dual  spectral 
bands  method,  it  is  required  to  show  one<to-one  correspondence  between  the  ratio  and 
the  temperature  of  the  object.  For  this  purpose,  the  derivative  of  the  ratio  of  band  lim¬ 
ited  photon  fluxes,  Sk,  in  (2.7),  with  respect  to  temperature  is  performed.  If  the  ratio  of 
photon  emittance  is  proved  to  be  either  strictly  increasing  or  strictly  decreasing 
function  with  respect  to  temperature,  i.e.,  either  d^JdT >  0  or  dStJdT  <  0  for  all  T,  then 
the  ratio  of  photon  emittances  and  the  temperature  exhibit  a  one-to-one  correspond¬ 
ence.  It  follows  that  this  approach  is  potentially  applicable  to  the  temperature  determi¬ 
nation  of  grey  bodies. 

To  simplify  the  integral  of  (2.2)  ,  the  variables  of  integration  in  photon  emittance 
Q{X„  Xj,T)  are  replaced  by  the  dimensionless  parameters  [Ref.  8:  pp.  1257] 

and 

a  =  ,  (3.2a) 

where  T,  is  an  arbitrary  reference  temperature.  By  prior  convention  [Ref.  8:  pp.  1257], 
it  is  chosen  to  be 


Tr  = 


he 


02b) 
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where  as  introduced  in  Chapter  IJ,  is  the  highest  wavelength  in  spectral  band  1.  Using 
(3.2b)  as  a  reference  temperature  in  (3.1),  the  spectral  band  limits  of  are  given  by; 


he 


(3.3) 


where  is  one  of  the  spectral  band  limits  [Ref.  8:  pp.  1257],  ARer  substitutions  into 
(2.5),  the  new  representation  for  the  photon  emittance  is  given  by 


QiX„  Xj,T)  -  2nck^T^c-^h-^r‘g{^)dil,, 


(3.4) 


where  [Ref.  8:  pp.  1257) 


«('/') 


expiik)  -  1  * 


The  ratio  of  photon  emittance  is  obtained  as 


(3.5) 


f^a 

J  * 
J  ^ 

•'l«3 


g{\tf)d\l/ 


gWd^ 


(3.6) 


and  the  derivative  of  the  ratio  of  photon  emittance  signals  (3.6)  with  respect  to  temper¬ 
ature  comes  out  to  be 
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where 


A^) 


[exp{\lf)  - 1]* 


(3.8) 


(The  derivation  of  (3.7)  is  in  Appendix  A.)  To  simpliiV  the  expression  of  (3.7),  the  fol« 
lowing  operator  is  defined  [Ref.  8:  pp.  1257J. 


Ia,6](l ,  jfj,  X3,  X4, «)  -  f  'a{tp)d\lt  f  -  f  'K'I'W  f  (3.9) 

•’00C3  •'eucj 


Then  the  ratio  becomes 


dT  ^  T 

) 

pax* 

2 

(3.10) 


Since  the  denominator  of  (3.10)  is  always  positive,  the  sign  of  dSlJdT  is  determined  by 
the  sign  of  the  operator  (f,g].  It  is  necessary  to  show  that  the  sign  of  [f,g],  which  is  the 
sign  of  (3.10),  does  not  change  with  «,  or  equivalently  with  temperature.  The  temper¬ 
ature  dependence  of  [f,g]  appears  only  in  oc,  which  can  be  seen  in  the  limits  of  integration. 
The  integrands  do  not  have  an  explicit  dependence  on  temperature. 

A.  NON-OVERLAPPING  CASE 

In  order  to  prove  that  the  sign  of  [f,g]  is  fixed  over  the  domain  of  a,  several  proper¬ 
ties  otA^)  required.  They  are  (Ref.  8:  pp.  1257] 

1. y(0)-^(0)-0 

2.  0  <  g{^l/)  <A^)  for  >  0 

3.  linv- J(^) "  li«V-o.?(^)  -  0 

4.  For  p  >  l,A^)  ~  PsW  has  exactly  one  root  on  the  interval  \l/  e  (0, 00). 
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The  discussion  to  follow  will  clarify  the  significance  of  these  requirements.  Grapliical 
representation  is  useful  to  clarify  these  properties  and  it  is  shown  in  Figure  2.  This  figure 
is  especially  useful  for  visualization  of  the  fourth  property.  For  values  of  /?  >  1,  it  is  ob* 
vious  thaty(^)  crosses  0gi\lt)  exactly  once.  A  more  formal  discussion  of  the  four  prop* 
erties  is  provided  in  the  next  paragraph. 

The  first  property  is  proven  by  using  THopital's  rule.  In  order  to  prove  the  second 
property,  define 


W- 


m 

gW  ' 


which  by  direct  substitution  from  (3.5)  and  (3.8)  is  given  by 


Pm 


1  -  cxp(  -  il/)  * 


(3.11a) 


(3.1 Ifc) 


It  then  can  be  shown  that  Pm  >  1  for  ^  >  0.  This  follows  from  the  derivative  of  P(i/'), 
which  is ; 


dPj^)  ^  1  -  cxp(  exp(  -  \j/) 

“  (I“Cxp(-i^))^ 


Afler  expanding  (3.12a),  it  is  clear  that : 

dPm  _  1 _ ^  exp(  -0)  ^  ^ 

d\lf  "  1  -  exp(  -  ^^r)  “  (1  _  cxp(  -  \li)f 


iZ.Ub) 


which  follows  from  the  observation  that  the  second  term  in  the  expansion  is  smaller  in 
magnitude.  Also,  by  direct  application  of  I'Hopital's  rule  on  (3.11b),  lin\^o+F(^)  •=  1. 
These  two  results  imply  that  P{\p)  >  1.  Thusy(iA)  >  Consequently,  after  noting  that 
bothyf^^'),  as  given  in  (3.8),  and  gm>  as  given  in  (3.5),  are  strictly  positive,  property  two 
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follows.  The  third  property  is  simply  proven  by  limit  evaluation.  For  the  fourth  prop¬ 
erty,  it  is  necessary  to  evaluate  -  0g{\lf)  given  by 


r 


icxpW-Xf  '^expW-l 


(3.13a) 


which  can  be  expressed  as 


Ji.vi  mvj  exp(rt-l  \^exp(|H-l 


(3.136) 


By  noting  (3.1  lb),  it  follows  that 


A^)  -  0gi'l') ' 


r 


—  LP(^^)~/?1 


(3.13c) 


Since  Iin\,o*fi('/')  “  1  0  is  constrained  to  be  greater  than  one,  the  quantity 

A^)  “  satisfies,  after  reference  to  (3.11a),  the  inequality 


A'lf)  -  0g{^)  <  0  (small 


(3.14a) 


for  suflicicntly  small  \j/.  Also  since  lin\^o«*^(l^)“oo  »  the  inequality  has  to  reverse  for 
large  enough  ilf.  In  other  words 


A'!')  -  Pgi^)  >  0  (large  \l;). 


(3.146) 


Since,  from  (3.12a),  P((/f)  is  a  strictly  increasing  function,  it  is  correct  to  conclude  the 
inequality  reversal  can  only  occur  once.  Property  4  is  therefore  proven. 

For  the  non-overlapping  case  represented  on  Figure  la,  there  are  only  two  sets  of 
relations  which  can  define  the  spectral  bands  (3.3): 


1  <  Jf2  ^  X3  <  X4, 


(3.15a) 
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and 


X3  <  X4  1  <  X3. 


(3.156) 


Relation  (3.1Sa)  is  considered  first.  Motivated  by  the  fourth  property  previously  stated, 
let 


«(«Jfa)  ’ 


(3.16) 


The  quantities  y(^)  and  P{oL)g{\lf)  must  cross  each  other  only  once  at  the  value  of  ^  de> 
noted  as  x,.  This  is  represented  in  Figure  3.  The  geometrical  aspects  of  figures  dictates 
the  following  relations  for  all  values  of  «: 


J. 


ouci 

MW  <  J  Pi»)MW 


(3.  Hu) 


and 


•'aotj  •'mj 


(3.176) 


Substitution  (3.17a)  and  (3.17b)  into  (3.9)  shows  the  operator  [f,g]  satisfies; 


I/‘^](l,Jf2.Jf3.<*'4.«)<  f  ^Pia)g{^w\*^gi^W-\  ’gi4>w\^Pi<!t)g{^W^0 

•'«  *’axj  *'a  •'otxj 


(3.18) 
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for  all  the  values  of  «  (3.2a).  According  to  the  brief  description  following  (3.10),  this 
proves  that  d0i,ldT<O  for  all  the  values  of  T  and  therefore  is  a  strictly  decreasing 
Ametion  of  temperature.  For  the  second  relation  of  non-overlapping  bands  (3.1Sb),  a 
similar  analysis  would  show  that  ddlJdT  >  0  and  Si,  is  a  strictly  increasing  function  of 
temperature. 

B.  PARTIALLY  OVERLAPPING  AND  COMMON  BAND  LIMIT  CASE 

In  this  section,  the  case  which  is  partially  overlapping  and  the  common  band  limit 
case  are  considered.  Two  sets  of  relations  are  considered  in  this  case  and  they  are 

1  ^  Jf3  <^2  <  Jf4  or  1  <  Jf3  <  JTj  S  Jf4  (3.19a) 

X3  s  1  <  jf4  <  jf2  or  jf3  <  1  <  jf4  jf2  (3. 19i>) 

For  this  discussion,  the  graphical  method  previously  introduced  is  useful  [Ref.  8:  pp, 
1258].  Figure  4  illustrates  the  first  set  of  the  relations  (3.19a)  for  partial  overlap.  With 
reference  to  Figure  4,  six  sets  of  integrals  are  defined  as  follows: 


PMj 

A]  “ 

•'a 

(3.20a,fc) 

f«l 

Am 

rout, 

(3.20c,</) 

•'flUCj 

Am 

Aj  - 

(3.20c,/) 

After  inserting  this  formalism  into  (3.9),  the  bracketed  operator  lf,g]  becomes 


IS 


^2*  -*'3»  •*4*  ®)  “  (^l  +  ^2)M2  "*■  ^3)  ”  (^i  +  '^2)('S>  +  S3). 


(3.2ia) 


After  expansion  and  recombination  of  (3.2 la),  it  follows 

(/•^Id.  Jf2.  Jf3.  Jf4. «)  -  M3(Si  +  Sj)  -  53(^1  +  ^2)]  +  (S,^j  -  ^jSj).  (3.216) 

This  can  be  retranslated  using  (3.20)  into  two  terms  which  follow 

[/',gKl.  Jf2.  Jf3.  *'4. «)  “  Jf2.  -*2*  *4. «)  +  t/i5l(l.  *3.  *3.  Xi,  a).  (3.2lc) 

Case  (3.19a)  will  now  be  applied  on  expression  (3  21c).  If  1  <X3<Xj<jc4,  then  both 
tetms  in  (3.21c)  are  of  the  form  1  <  jT}  «  jtj  <  jr4  which  satisfies  the  nonovertap  condition 
(3.15a).  Since,  as  indicated  in  the  previous  section,  the  particular  nonoverlap  condition 
(3.15a)  generates  negative  values  for  (f,gl,  the  expression  (3.21c)  must  be  negative  for  the 
stated  inequality. 

For  common  band  limit  case,  which  is  1  -•  jc,  <  x,  <  X4  or  1  <  Xj  <  Xj » jr4 ,  one  of  the 
terms  in  (3.2Ic)  becomes  zero.  However,  the  other  term  remains  negative,  and  therefore 
(f,g]  remains  negative.  This  discussion  implies  that  the  ratio  01^  is  a  strictly  decreasing 
function  of  temperature.  Similarly,  for  the  second  relation  of  this  case  (3.19b),  the  ratio 
is  a  strictly  increasing  function  of  temperature. 

C.  TOTALLY  CONTAINED  SPECTRAL-BAND  CASE 

This  section  discusses  case  3,  where  one  interval  is  entirely  contained  within  the 
other.  By  the  use  of  asymptotic  approximations  [Ref.  8:  pp.  1258],  it  is  possible  to  show 
that  the  ratio  Si,  may  have  the  same  value  for  two  or  more  temperatures.  Therefore,  the 
measurement  to  determine  the  temperature  is  not  unique.  Two  examples,  which  are  the 
case  where  a  is  large  and  the  case  where  a  is  small,  are  discussed. 

As  previously  discussed,  the  sign  of  dSiJdT is  determined  by  the  sign  of  [f,g].  There 
are  two  sets  of  relations  which  are  1  <  x,  <  x*  <  x,  and  Xi  <  1  <  x,  <  X4.  The  first  set  of 
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relations  is  arbitrarily  chosen  for  discussion  since  the  approach  is  similar  for  the  other 
relation.  For  large  or  small  value  of  a  asymptotic  methods  are  useful  to  determine  the 
sign  of  [f,g].  (The  derivation  using  asymptotic  approximation  is  in  Appendix  B.)  For  the 
large  value  of  a  (i.e.,  small  temperature  ),  (f,g]  is  approximated  as 

-  X3)  expC  -  «(X3  +  1)].  (3.22) 

Since  jtj  is  greater  than  one,  [f,g|  becomes  negative.  For  small  values  of  a  (i.e.,  large 
temperature  ),  [f,g]  is  approximated  as: 

-fj  ~  -  4)1  (3.23) 

Given  and  jc'4,  it  is  always  possible  to  find  a  sufficiently  large  x^  such  that  (3.23)  is 
positive.  Combining  the  tw'o  results  for  small  a  and  large  a,  the  sign  of  d^JdT  is  not 
fixed  over  a  (i.e.,  temperature  ).  This  means  that  the  ratio  does  not  have  a  one-to-one 
correspondence  with  temperature.  Therefore,  for  the  case  when  one  spectral  band  is 
completely  contained  within  the  other,  the  temperature  determination  method  based  on 
the  dual  spectral  band  measurement  of  photon  emittance  is  not,  in  general,  valid. 

D.  RATIO  CALCULATION  FOR  PHOTON  EMITTANCE 

In  order  to  generate  the  output  ratio  calculation,  the  numerical  method  for  inte¬ 
gration  known  as  Simpson's  rule  (Ref.  12:  pp.  95]  was  used.  The  flow  chart  for  this 
calculation  is  shown  in  Appendbc  C.  The  examples  were  demonstrated  using  the  condi¬ 
tions  in  [Ref  8:  pp.  1258-1259).  One  case  has  a  band  combination  of  (2.5  nm,  5.5  /mi] 
and  (3.9  nm,  4.1  //m]  shown  in  Figure  5.  The  second  case  has  a  band  combination  of 
(10.0  nm,  14.4  /im]  and  (11.8  fim,  12.2  /im]  shown  in  Figure  6.  In  these  figures,  the  nor¬ 
malized  ratio  denoted  as  was  used  to  enhance  the  results.  Normalization  is  done  by 
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(3.24) 


~  ^mln 

- sr- 


mm 


-M'here  is  the  maximum  value  of  St,  and  is  the  minimum  value  of  dt,  Both  cases 

are  slightly  different  from  the  results  in  [Ref.  8:  pp  1259],  but  the  general  features  are  the 
same.  One  ratio  value  may  correspond  to  two  or  more  temperature  values  as  predicted 
in  Section  C.  These  figures  indicate  that  the  minima  of  the  ratio  occur  in  different  po¬ 
sitions  depending  on  the  bands  selected. 

E.  DISCUSSION  OF  COMPUTER  RESULTS 

Computer  simulations  were  performed  for  the  three  cases  discussed  in  this  chapter. 
The  first  example  is  the  nonoverlapping  case  and  the  band  pair  is  chosen  to  be  [3.0  fim, 
4.0  nm]  and  [4.0  ^un,  5.0  nm\.  Figure  7  shows  this  example  and  it  exhibits  the  decreasing 
feature  as  expected  in  Section  A.  The  second  example  is  the  partially  overlapping  case 
and  the  band  pair  is  chosen  to  be  [3.0  fim,  4.5  ^ni]  and  [3.5  (im,  5.0  /itn).  Figure  8  shows 
this  case  and  the  results  agree  with  the  discussion  in  Section  B.  The  third  example  is  the 
case  of  totally  contained  band  within  the  other  and  the  band  pair  is  [3.0  ixm,  5.0  mhi]  and 
[3.5  nm,  4.5  Mtn].  Figure  9  represents  this  case  and  the  results  show  that  one-to-one 
correspondence  does  not  exist  in  this  case.  It  can  be  seen  from  Figure  9  that  a  ratio  value 
of  2.0  corresponds  to  both  460  degrees  Kelvin  and  1400  degrees. 
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Figure  3.  Nonoverlapping  case.:  The  constant  j?(a)  is  chosen  such  that 

0{<l)g(<tXi)  -  flojf,). 


Figure  4.  Partially  overlapping  case  or  common  band  limit  case:  Each  region  is 
defined  in  the  text. 
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Figure  5.  Example  of  totally  contained  case  (normalized  ratio):  Bands  are  [2.5 
Mm,  5.5  Mm]  and  [3.9  Mm,  4.1  Mm]. 
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Figure  6.  Example  of  totally  contained  case  (normalized  ratio):  Bands  are  [10 
Aim,  14.4  Aim]  and  [11.8  Aim,  12.2  Aim]. 


f  igure  7.  Example  of  the  non-overlapping  cose:  Bands  are  [4.0  nm,  5.0  nm]  and 
[3.0  tim,  4.0  /«n]. 
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Figure  8.  Example  of  the  partially  overlapping  case:  Bands  are  [3.5  /xm,  5.0  /xm] 
and  [3.0  ^m,  4.5  /xm]. 
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Figure  9.  Example  of  the  totally  contained  case:  Bands  are  [3.0  nm,  5.0  /im]  and 
[3.5  Atm,  4.5  /xm]. 
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IV.  A  GENERAL  MODEL  FOR  TEMPERATURE  DETERMINATION  VIA 
DUAL  SPECTRAL  BAND  MEASUREMENTS 

In  Chapter  II,  ideal  photon  emittance  temperature  determination  was  discussed. 
Since  this  discussion  was  dealing  with  the  ratio  without  considering  measurement  of  the 
photon  emittance  as  a  system,  it  was  not  always  practical.  Therefore,  some  other  factors 
should  be  taken  into  account  to  realize  this  method.  They  are  emissivity  of  the  object, 
atmospheric  transmittance,  detector  responsivity,  and  so  on.  In  the  previous  discussion, 
those  parameters  were  assumed  to  be  spectrally  constant  or  unity.  Strictly  speaking, 
they  depend  on  the  wavelength.  Also  after  considering  these  factors  as  part  of  the  sys¬ 
tem,  the  resultant  ratio  becomes  the  ratio  of  the  output  voltages  from  the  detectors  of 
two  bands  instead  of  the  ratio  of  photon  emittance.  A  brief  discussion  of  these  param¬ 
eters  follow. 

A.  EMISSIVITY 

The  emissivity  dictates  how  much  thermal  energy  is  radiated  from  an  object.  Its 
value  lies  between  zero  and  unity  depending  upon  the  material  characteristics  of  thermal 
source.  It  is  denoted  as  e(>l).  Although  emissivity  for  some  materials  has  temperature 
dependency,  it  is  assumed  to  be  independent  of  temperature  to  make  discussion  simpli¬ 
fied.  Materials  are  divided  into  three  types  based  on  the  value  of  emissivity.  They  are 

(Ref.  3:  pp.  40] 

1.  A  black  body,  c(.A)  =  1. 

2.  A  gray  body,  c(ii) »  constant  <  1. 

3.  A  selective  absorber,  c(/l)  varies  with  wavelength. 

The  spectrally  constant  emissivity  cases  were  covered  in  the  previous  Chapters  II  and 
III.  Therefore,  wavelength-dependent  emissivity  is  the  primary  focus  of  discussion  in  this 


chapter.  In  addition,  since  emissivity  depends  on  the  types  of  materials,  it  is  not  possible 
to  get  reliable  analytical  expressions  for  all  selective  absorbers.  Usually  emissivity  values 
are  obtained  by  measurement. 

B.  ATMOSPHERIC  TRANSMITTANCE 

The  atmosphere  is  known  to  behave  as  a  wavelength-selective  absorber.  This  im¬ 
plies  that  the  transmission  of  signals  through  the  atmosphere  becomes  wavelength  de¬ 
pendent.  Specifically,  the  atmosphere  affects  the  transmission  of  infrared  radiation.  The 
main  compositions  of  the  atmosphere  are  nitrogen,  carbon  dioxide,  oxygen,  and  water 
vapor.  Among  these,  molecules  which  have  nonsymmetric  molecular  structure  are  pri¬ 
mary  responsible  for  absorption.  Thus  carbon  dioxide  and  water  vapor  are  the  dominant 
sources  for  absorption.  In  addition  some  other  small  constituents  of  the  atmosphere, 
which  have  non-symmetric  structure,  also  contribute  absorption. 

Atmospheric  scattering  is  the  other  important  phenomenon  which  affects  atmo¬ 
spheric  transmission.  The  strongest  atmospheric  scattering  occurs  in  the  case  when  the 
size  of  aerosols  and  the  wavelength  of  the  signal  are  almost  the  same.  In  this  particular 
case  the  scattering  coefficient  depends  highly  on  a  wavelength  of  the  signal.  The  Mie 
theory  describes  this  phenomenon. 

The  atmospheric  transmittance  is  determined  by  atmospheric  absorption  coefficient 
and  atmospheric  scattering  coefficient.  It  is  given  by 

T(.l,j:)-cxp[-(a(/l)-l-y(.l))x]  (4.1) 

where  a(ii)  is  atmospheric  absorption  coeflicient,  y(^)  is  atmospheric  scattering  coeffi¬ 
cient,  and  X  is  the  distance  between  the  object  and  the  detector.  An  example 
transmittance  data  by  the  NPS  LOWTRAN6  program  [Ref.  13]  at  20  meters  in  height, 
0  degree  elevation  for  1  km  and  5  km  distances  under  1962  U.S.  standard  atmospheric 
conditions  arc  shown  in  Figure  10.  The  importance  here  is  that  the  transmittance  de- 


pends  not  only  the  wavelength  but  also  distance  and  atmospheric  conditions.  Moreover, 
atmospheric  conditions  vary  with  height  and  thus  the  transmission  depends  on  the  path 
between  the  object  and  the  detector.  Unless  the  path  and  atmospheric  conditions  arc 
known,  it  is  impossible  to  apply  the  ideal  ratio  calculation.  Therefore,  the  transmission 
path  and  atmospheric  conditions  have  to  be  provided  in  order  to  obtain  the  atmospheric 
transmittance.  A  recent  report  (Ref.  1 1:  pp.  41]  on  transmittance  calculation  for  the  at¬ 
mosphere  has  concluded  that  refraction  effects  can  essentially  be  ignored. 

C.  DETECTOR  RESPONSIVITY 

The  detector  responsivity  is  generally  obtained  by  dividing  the  measured  output 
signal  voltage  by  the  optical  power  falling  on  the  detector.  And  it  is  defined  as 


R  s 


Kms 

HA, 


(4.2) 


where  R  is  the  responsivity  [V/W],  V,„,  is  the  RMS  output  voltage,  H  is  the  RMS 
irradiance  [IK/cm’j  at  the  detector  location,  and  is  the  detector  area  [cnt^].  However, 
this  quantity  is  not  the  appropriate  parameter  for  the  development  of  the  photon 
emittance  case.  A  new  proposed  responsivity  R'  is  defined  by  modifying  (4.2)  in  order 
to  adjust  it  for  the  photon  emittance  calculation.  It  is  given  by  dividing  the  measured 
output  voltage  signal  by  the  number  of  photons  into  the  detector.  The  proposed  R'  is 
defined  as 


R' 


(4.3) 


where  R'  is  the  photon  responsivity  [V •  photonir' •  s],  and  Q  is  the  photon  emittance 
[photons  •  r' •  cm-^.  It  follow's  from  definitions  (4.2)  and  (4.3)  that  the  photon 
responsivity  R'  is  simply  obtained  by  multiplying  responsivity  R  by  hcj?., 
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(4.4) 


as  shown  above.  For  photon  detectors  sensitive  to  the  number  of  photons,  R'  should 
be  fairly  constant  out  to  the  cut  off  wavelength. 

Unfortunately,  detector  responsivity  (4.2)  is  not  usually  provided  among  detector 
specifications,  but  the  specific  detectivity  D*  is  given  instead.  The  D*  is  the  normalized 
quantity  of  detectivity  defined  by  [Ref.  3:  pp.  270]: 

D*  -  (4.5a) 


where  B  is  the  electrical  system  bandviidth  in  heltz,  A4  is  the  area  of  the  detector,  and 
D  is  the  detector  detectivity.  This  quantity  is  oAen  expressed  in  terms  of  the  reciprocal 
of  the  noise  equivalent  power  (NEP).  Therefore,  the  D*  is  also  given  by  [Ref.  3:  pp. 
270] 

(4.5« 


The  noise  equivalent  power  is  the  radiant  flux  required  to  obtain  an  output  signal  equal 
to  the  noise  of  the  detector.  It  is  given  by 


NEP 


HAg 

Vn 


(4.6) 


where  V,  is  RMS  signal  voltage  and  is  RMS  noise  voltage.  From  (4.6),  the  specific 
detectivity  becomes 


D* 


HAa 


(4.7) 
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The  photon  responsivity  is  obtained  by  combining  (4.3),  (4.4),  and  (4.7),  and  it  is 


R' 


D*Vn  he 
yfAj  ^  ‘ 


(4.8) 


This  equation  shows  the  specific  deteaivity  is  closely  related  to  the  photon  responsivity. 
From  the  expected  wavelength  dependencies  in  (4.8),  it  follows  that  the  quantity  D* 
should  increase  linearly  with  wavelength  out  to  the  cutoff  wavelength  (Ref,  3:  pp.  294]. 
Since  photon  detectors  respond  to  the  number  of  photons  absorbed  in  the  detector  up 
to  cutoff  wavelength  (usually  dictated  by  the  band'gap),  the  modified  responsivity  of 
photon  detectors  is  ideally  independent  of  wavelength. 

D,  COMPUTER  SIMULATION  OF  TEMPERATURE  DETERMINATION 

Based  on  above  discussion,  the  ratio  of  photon  emittance  for  black  body  is  replaced 
by  the  ratio  of  the  output  voltages  at  the  detectors,  and  it  is  given  by 


J%A)T(A,x)R',(;)g(^)rf^ 

-  m 


where  c(A)  is  the  emissivity,  t(A,  x)  is  atmospheric  transmittance,  and  R'(A)  is  the  photon 
detector  responsivity.  The  subscript  on  the  photon  detector  responsivity  denotes  the 
spectral  band  used  in  the  measurement.  In  this  equation,  assuming  that  the  wavelength 
parameters  are  spectrally  constant  and  two  identical  detectors  are  used,  the  ratio  be¬ 
comes  identical  to  the  ratio  of  photon  emittance  for  black  body  (3.6).  In  order  to  gen¬ 
erate  the  output  ratio  calculation  under  more  general  conditions,  the  numerical  method 
for  integration  known  as  Simpson's  rule  [Ref.  12:  pp.  95]  is  used.  The  algorithm  for  this 
calculation  is  shown  in  Appendix  D. 
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1.  Discussion  of  the  computer  simulations 

As  described  in  this  chapter,  the  actual  ratio  for  the  photon  detector  is  the  ratio 
of  the  output  voltages.  In  the  rest  of  this  section,  the  simulation  for  the  ratio  calculation 
is  discussed.  Although  the  photon  detector  responsivity  has  wavelength  dependence,  as 
discussed  in  Section  C,  it  is  not  a  very  sensitive  factor  up  to  the  detector  cut-off. 
Therefore,  this  is  ignored  in  this  discussion,  and  emissivity  and  atmospheric 
transmittance  are  considered.  Since  emissivity  depends  on  the  material  of  the  thermal 
source,  many  kinds  of  emissivity  values  are  possible.  For  this  simulation,  an  arbitrary 
created  emissivity  curve  shown  in  Figure  11  was  used.  For  atmospheric  transmittance, 
the  atmospheric  conditions  of  1962  U.S.  standard  atmosphere  were  selected.  The  data, 
shown  in  Figure  10,  are  obtained  from  NFS  LOWTRAN6  program.  For  the  purpose  of 
comparison  with  the  results  for  ideal  case,  the  same  band  pairs  are  selected  as  ones  in 
Chapter  III.  The  result  for  overlapping  case  is  shown  in  Figure  12,  for  partially  over¬ 
lapping  case  in  Figure  13,  and  totally  contained  case  in  Figure  14.  The  curves  generated 
in  these  figures  correspond  to  1  km  and  5  km.  With  the  exception  of  Figure  14,  there 
is  a  one-to-one  correspondence  between  measured  ratio  and  temperature  for  each  range. 
In  Figure  14,  there  is  a  slight  positive  shape  after  T*1300  K  which  violates  the 
uniqueness  condition.  Although  the  values  arc  shifted  slightly,  the  shapes  of  the  curves 
are  essentially  the  same  as  ones  for  ideal  ratio  case  shown  in  Chapter  III.  The  general 
features  of  each  case  are  the  same.  However,  the  results  demonstrate  that  the  curves 
depend  on  distance.  This  is  quite  inconvenient.  The  application  of  two-band  temperature 
determination  is  then  impossible  without  range  information.  Also,  to  lesser  extent,  the 
application  is  limited  since  the  emissivity  is  material  dependent.  If  the  source  for  thermal 
radiation  is  not  a  gray  body,  then  the  emissivity  of  the  source  needs  to  be  known  in  ad¬ 
vance. 
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Figure  10.  Atmospheric  transmittance  for  1962  U.S.  standard  atmosphere:  The 
path  lengths  are  taken  to  be  1  km  and  5  km  in  the  horizontal  direction. 
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Figure  1 1.  Example  of  the  emissivity  (arbitrary  chosen):  This  data  does  not  cor¬ 
respond  to  the  specific  material. 
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Figure  12.  The  ratio  of  output  voltages  for  nonoverlapping  case:  Bands  are  (5.0 
nm,  4.0  fim)  and  (4.0  nm,  3.0  ^m).  The  distances  are  1  km  and  5  km. 


35 


37 


V.  THERMAL  RANGING  TECHNIQUE 


From  the  previous  chapter,  the  ratio  of  output  voltages  from  dual  spectral  bands  is 
not  the  same  as  the  ratio  of  black  body  photon  emittances  under  realistic  and  nonideal 
conditions.  To  solve  and  apply  the  method  under  these  conditions  can  be  quite  difficult. 
On  the  other  hand,  one  signiiicant  idea  is  hidden  in  this  problem. 

The  concept  presented  in  this  chapter  is  created  from  the  difficulty,  in  the  temper¬ 
ature  determination  method,  produced  by  nonconstant  atmospheric  transmittance. 
Among  assumptions  discussed  in  Chapter  IV,  atmospheric  transmittance  breaks  the  va¬ 
lidity  of  temperature  determination  procedure  because  it  depends  on  wavelength,  the 
distance  between  an  object  and  the  detector,  and  atmospheric  conditions.  Although 
emissisivity  and  detector  responsivity  have  to  be  taken  into  account,  under  some  rea¬ 
sonable  assumptions  they  can  be  neglected.  It  will  be  shown  that  the  dependency  of  at¬ 
mospheric  attenuation  on  the  distance  makes  the  concept  of  thermal  ranging  feasible. 

The  use  of  infrared  thermal  signals  has  the  advantage  that  active  interrogation  is  not 
required  since  the  object  radiates  energy  or  photons.  This  passive  property  of  infrared 
detection  can  be  significantly  important  in  military  situations.  However,  the  passive  de¬ 
tection  of  the  target  does  not  automatically  provide  range  data  and  this  is  a  disadvantage 
of  thermal  infrared  detection  systems.  One  interactive  (nonautomatic)  passive  system 
involves  the  use  of  a  trained  human  observer.  The  training  allows  the  observer  to  rec¬ 
ognize  a  target  and  by  noting  the  size  of  the  target  on  the  screen  make  an  experienced 
estimate  of  the  range.  Currently  ranging  is  done  by  active  sensor  systems  such  as  radar 
or  laser  range  finder.  These  active  systems  are  capable  of  being  sensed  by  the  target. 
Although  active  ranging  systems  tend  to  provides  precise  distance  information,  they  ex- 
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pose  the  instrument  generating  the  signal.  Therefore,  the  development  of  a  passive 
ranging  system  would  provide  significant  military  advantage. 

The  Electronic  Support  Measure  (ESM)  equipment  in  electronic  warfare  systems  is 
the  passive  sensor  in  the  microwave  region.  This  gives  target  identification,  direction  of 
the  target  and  also  a  very  rough  range  information.  For  this  method,  the  power  trans¬ 
mitted  by  a  radar  of  a  target,  the  antenna  gain,  and  attenuation  factor  must  be  obtained 
in  advance.  It  seems  that  the  same  concept  can  be  applied  to  infrared  systems  only  if 
the  temperature  of  a  target,  surface  area,  emissivity,  and  atmospheric  transmittance  are 
known.  Therefore,  this  method  is  not  practical. 

The  follo^^'ing  discussion  proposes  a  new  concept  of  passive  ranging  by  infrared 
systems.  It  will  be  shown  that  this  concept  is  a  feasible  solution  to  the  problem  of  target 
ranging  using  only  automatic  passive  methods.  The  following  discussion  consists  of 
three  steps: 

1.  Step  1:  Atmospheric  attenuation  coefficient  is  assumed  to  be  constant.  Temper¬ 
ature  of  the  target  is  known. 

2.  Step  2:  Atmospheric  attenuation  coefficient  is  assumed  to  be  constant.  Temper¬ 
ature  of  the  target  is  unknown. 

3.  Step  3:  Atmospheric  attenuation  coefficient  is  not  constant.  Temperature  of  the 
target  is  unknown. 

The  concept  of  ranging  shown  here  is  based  on  the  distance  dependency  of  atmospheric 
transmittance.  Step  1  and  step  2  are  intermediate  procedures  and  they  do  not  allow 
changing  atmospheric  conditions  due  to,  for  example,  an  elevated  transmission  path. 
In  step  3,  these  prohibitions  are  rescinded  and  the  ranging  model  has  greater  applica¬ 
bility. 

A.  RANGING  FOR  KNOWN  TEMPERATURE  CONDITION  WITH 
HOMOGENEOUS  ATMOSPHERIC  CONDITIONS 

In  the  analysis  of  this  section  the  emissivity  of  the  target  and  photon  detector 
responsivity  are  assumed  to  be  spectrally  constant.  In  addition,  it  is  assumed  that  at- 
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mospheric  'londitions  are  not  changing  along  a  line^of-sight  trajectory  between  target 
and  passive  sensor  system.  Therefore,  atmospheric  attenuation  coefficients  will  not  de* 
pend  on  range  but  still  exhibit  wavelength  dependence.  By  appropriate  band  selection, 
i.e.,  from  a  flat  portion  of  the  transmission  curve  [LOWTRAN],  the  wavelength  de¬ 
pendence  of  the  transmission  can  be  approximately  neglected  with  in  each  of  the  selected 
bands.  Therefore,  atmospheric  transmittance  for  the  i-th  band  is  given  by 

TXx)-exp(-<ypf)  (5.1) 

where  a,  is  the  atmospheric  attenuation  coefficient  within  the  /-th  band,  and  x  is  the 
transmission  path  length.  Applying  the  assumptions  stated  to  the  ratio  (4.9)  leads  to 

T2(jf) 

Afler  substitution  of  (5.1),  the  ratio  can  be  written  as 

5P*'-exp{-(a,-iTj)x]5?^  (5.26) 


where  dl**  is  the  ratio  of  the  measured  output  voltage  at  the  detector,  is  the  ideal  ratio 
of  photon  emittance,  and  a,  are  the  atmospheric  attenuation  coefficient  for  the  measured 
bands.  In  (5.2b),  the  ratio  of  the  output  voltages  is  obtained  by  direct  measurement.  The 
atmospheric  attenuation  coefficient  is  calculated  from  the  atmospheric  conditions,  and 
the  ratio  of  photon  emittance  is  determined  by  the  temperature  given.  Therefore,  the 
transmission  path  length  is  the  only  unknown  parameter.  From  (5.2b),  the  transmission 
path  distance  is  obtained  as 
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£ 


(5.3) 
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The  use  of  (5.3)  makes  passive  thermal  ranging  feasible.  However,  it  is  applicable  to  very 
limited  conditions.  In  the  battle  field,  the  targets  are  usually  unknoN^m,  and  consequently 
the  temperature  of  the  target  is  hard  to  obtain.  This  means  that  the  system  should  be 
capable  of  determining  both  range  and  temperature  to  be  practical.  The  new  proposed 
concept  is  described  in  the  following  sections. 


B.  RANGING  FOR  UNKNOWN  TEMPERATURE  WITH  HOMOGENEOUS 
ATMOSPHERIC  CONDITIONS 

In  this  section  a  concept  to  determine  the  range  and  the  temperature  of  the  target 
is  discussed.  As  in  Section  A,  the  atmospheric  conditions  are  assumed  to  be  homogene¬ 
ous  in  space,  Mathematically,  two  nondegenerate  equations  can  be  used  to  solve  for  two 
unknowns.  In  this  discussion,  the  two  unknowns  are  range  and  temperature  of  a  target. 
To  get  two  equations,  tw'o  distinct  sets  of  ratio  calculations  are  required.  For  that  pur¬ 
pose,  no  less  than  three  sets  of  band  measurements  ( i.e.,  band  1,  band  2,  and  band  3  ) 
are  required  to  produce  two  distinct  equations  such  as  (5.3).  All  the  assumptions  de¬ 
scribed  in  the  previous  section  are  also  applied  in  this  section  except  that  the  temper¬ 
ature  of  a  target  is  unknown.  The  mathematical  derivation  to  get  the  range  equations 
is  exactly  the  same  as  in  the  previous  section.  The  two  range  equations  for  three  sets  of 
band  measurements  are 
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(5.5) 


where  the  subscripts  serve  as  an  index  to  identify  spectral  bands.  Combining  (5.4)  and 
(5.5),  one  equation  to  solve  the  temperature  is  obtained  as 
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(5.6) 


1 


u 


1  ^23 


However,  tlris  equation  is  impossible  to  solve  analytically.  A  numerical  method  can  be 
used  to  solve  this  equation.  After  getting  the  temperature,  the  range  is  easily  obtained 
by  using  (5.4)  or  (5.5).  Since  the  assumption  that  the  atmospheric  conditions  are  ho> 
mogeneous  in  space  is  only  approximately  valid  for  very  short  ranges,  this  method  is  not 
practical.  The  following  section  solves  the  difTiculty  of  the  determination  of  the  range 
without  assuming  homogeneous  atmospheric  conditions. 


C.  RANGING  FOR  UNKNOWN  TEMPERATURE  WITH  INHOMOGENEOUS 
ATMOSPHERIC  CONDITIONS 

For  the  following  discussion,  several  equations  are  recalled  for  convenience,  For 
four  bands,  the  ideal  ratios  calculated  from  photon  emittance  are  defmed  as; 


,/  ^2»  ^ 


and 


where  the  total  photon  emittance  is  given  by 


(5.7a) 


iS.lb) 


(5.8) 


and  spectral  photon  emittance  Qt  can  be  expressed  in  terms  of  the  gray  body  distribution 
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(5.9) 


Qa 


2ntc 


r{exp(/icM70-lJ 


With  two  band  measurements,  the  temperature  can  be  determined.  With  three  band 
measurements,  it  should  be  possible  to  determine  both  temperature  and  range.  As  pre« 
viously  discussed,  range  information  is  available  as  long  as  the  effective  atmospheric  at> 
tenuations  in  each  band  are  distinct.  The  problem  posed  by  obtaining  range  information 
for  a  thermal  source  at  the  same  elevation  as  the  receiver  (homogeneous  atmospheric 
conditions)  is  not  as  complex  as  that  associated  with  arbitrary  elevation.  Elevation 
changes  in  the  trajectory  of  the  radiation  usually  require  consideration  due  to  the  cur< 
vature  of  the  earth.  In  the  general  problem  to  be  solved,  both  the  receiver  and  the  ther¬ 
mal  source  take  arbitrary  locations  in  the  plane  defined  by  the  vertical  line  passing 
through  the  receiver  and  thermal  source.  This  is  represented  in  Figure  IS. 

In  the  algorithm  to  be  discussed,  three  arbitrary  but  distinct  radiation  bands  are  se¬ 
lected.  Initially  disregarding  the  effects  due  to  the  atmospheric  attenuation,  calculated 
ideal  ratios  of  signals  are  given  by 


Ideal  signal  Band  1 
Ideal  signal  Band  2 


(5.10a) 


Ideal  signal  Band  3 
Ideal  signal  Band  1 


(5.106) 


Figure  16  shows  the  temperature  dependence  of  these  curves  for  bands  defined  by; 

,12)  ~  (i2.0  fim,  10.8  urn) 


(^3.  -^-a)  **  (iO-8  m,  9.6  urn) 
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(^5«  -*6)  •  (9-6  MW,  8.2  urn) 


The  atmospheric  attenuation  in  the  form  of  transmittance  is  provided  from  the  NPS 
LOWTRAN6  program  [Ref.  13]  using  data  corresponding  to  a  vertical  path.  Atten¬ 
uation  coefficients  are  derived  from 

In  [nepersim]  (5. 1 1) 

where  i- 1,2,3  indicate  the  wavelength  band  and^‘  is  the  index  of  the  specific  layer.  The 
quantity  Tj  is  the  length  of  the  trajectory  in  the^-th  atmospheric  layer.  See  Figure  17. 
Each  fj  has  the  same  length  if  the  line-of-sight  angle  is  90  degree.  The  attenuation  coef¬ 
ficient  is  obtained  by  taking  the  90  degree  angle.  The  standard  band  average  for 
is  given  by 


T(/» 


(5.12) 


With  reference  to  Figure  18,  the  effective  approximate  transmission  factor  t„  for  the  /-th 
band,  is  given  by; 


^^nut 

T,  -  exp(  -  (5. 1 3) 

where  is  the  number  of  layers  between  source  and  receiver,  and  m  is  the  layer  num¬ 
ber  of  the  receiver.  Here,  r*  is  the  length  of  the  trajectory  in  the  layer  k  +  nt  -  I .  The 
approximation  involved  in  (5.12)  will  be  discussed  at  the  end  of  this  section.  Expression 
(5.13)  can  be  recast  into  the  form; 
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t,-exp(-a,R) 


(5.14) 


where 


kmi _ 

R 


(5.15) 


and  when  the  range  R  must  satisfy: 


(5.16) 


With  the  exception  of  the  first  and  final  layer,  the  value  for  each  r,  is  obtained  from  the 
rule 


n- 


(5.17) 


as  shown  in  Figure  19.  Each  R,  is  calculated  from  the  Cosine  Law  relation  derived  from 
the  geometry  of  Figure  20 

{H,  +  Re?  -  R}  +  {h^  +  Re?  -  IR^Jr^  +  Re)  cos(6  +  -f)  (5. 1 8a) 

or  equivalently 

{Hi  +  Re?  -  Rf  +  (/Jrf  +  Re?  +  2R,{hd  +  Re)  sin  6,  (5.1  U) 
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where  B  is  the  line  of  sight  angle  between  receiver  and  source  measured  with  respect  to 
the  horizontal.  R,  can  be  obtained  from  (5.18)  using  the  quadratic  formula  with  the 
correctly  chosen  '  +  '  sign.  /?,  is  also  obtained  from  (5.18)  after  taking 

+  (5.19) 

where  the  int[J  results  in  the  integer  part.  This  follows  from  examination  of  the  geom¬ 
etry  in  Figure  15.  For  R^  values  for  //*  follow 

//*«//,  +  (*-  l)Ar  (5.20) 

which  are  used  in  (5.18)  to  calculate  subsequent  For  clarity,  it  is  noted  that; 

/?,  -  r,  (5.21fl) 

»■  Tj  +  Ti  "f....  -f-  r*  (5.21ib) 

until  k  -  as  defined  by  (5.16).  Equations  (5.17)  combined  with  (5.20)  are  solved  it¬ 
eratively  to  generate  the  values  for  r,  except  for  r^,,.  Lastly,  is  calculated  from: 

-Rk^^-y  (5*22) 

Once  r,  are  determined  for  all  /  up  to  >  the  effective  transmission  for  each  band  (5.14) 
can  be  calculated  from  (5.15)  and  (5.16).  The  measured  ratios,  i.e.,  including  the  effects 
of  the  atmosphere,  are  related  tc  ideal  ratios  (5.7)  according  to: 

(5.J30) 

and 
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2l 

4^31  ■•4^31  T,  » 


(5.236) 


in  agreement  with  (5.2a).  AAer  substitution  from  (5.14),  (5.23)  can  be  written  as 

^^-^{2exp[~(<T,-<T2)-Rl  (5-24«) 


and 


^3\  “  ^31  I  -  (®3  “ 


(5.246) 


where  R  is  the  distance  between  thermal  source  and  receiver.  Figure  21  shows  the 

measured  ratios  calculated  from  (5.24)  for  conditions; 

•  O',  ■■  0.0001494  [nepersim] 

•  Oj  «■  0.0001027  [nepersim] 

•  03  *  0.0004622  [nepersim] 

•  R "  100  [m] 

Solving  for  R  from  (5.24),  it  follows  that 


R^\T) 


1 

O,  —  O2 


In 


{5.25a) 


and 


R^\T) 


_ 


(5.256) 


where  the  subscripts  are  employed  to  denote  spectral  bands.  In  the  algorithm,  it  will  be 
assumed  that  the  correct  measured  ratios,  consistent  with  (5.24),  arc  provided  by  exper¬ 
imental  testing.  Ideal  ratio  curves  can  always  be  calculated  according  to  the  methods 
(5.10)  previously  discussed.  Constants  a,,  Oj,  can  be  approximately  calculated  from  the 
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NPS  L0WTRAN6  data.  By  checking  for  the  crossing  point  between  the  two  curves 
(5.25),  both  the  temperature,  T,  and  the  range,  R,  are  uniquely  determined  as  shown  in 
Figure  22. 

The  algorithm  sketched  out  is  not  yet  complete.  The  values  of  a,,  (t,,  and  9,  are  de> 
pendent  on  the  range,  R,  since  the  weights  for  evaluating  the  coefficients  a,  (5.15)  depend 
on  all  the  values  for  r,.  In  order  to  calculate  the  coefficients  a,  (5.15),  an  initial  guess  is 
made  for  a  range  candidate,  /{«,,  within  limits  dictated  by  the  domain  of  LOWTRAN 
data  sampled.  This  value  for  Rg^  is  used  in  the  calculation  of  oi,  9},  and  93.  The  search 
for  a  numerical  crossing  point  for  curves  (5.25)  produces  a  new  value  for  R  which  can 
be  denoted  as  Rgi.  The  value  for  Rci  generated  by  this  first  iteration  is  expected  to  be 
closer  to  the  actual  range  but  with  an  unacceptable  level  of  error.  The  process  is  re¬ 
peated  by  iterating  with  the  R^.  The  next  estimate  generated  by  the  algorithm,  is 
tested  against  the  previous  one,  Rgy  In  general,  the  procedure  is  stopped  when 


<6. 


(5.26) 


where  d  is  the  predetermined  convergence  parameter  and  the  index  m  specifies  the 
number  of  iteration  cycles.  In  Appendix  E,  the  algorithm  for  the  process  is  presented. 

For  the  purpose  of  checking  the  algorithm,  the  range  candidate,  R(„,  is  selected  as 
shown  in  Figure  23.  Actual  ranges  are  chosen  from  10  meters  to  700  meters  and  Rga’s 
are  700  meters  for  up  to  350  meters  of  actual  ranges  and  10  meters  for  longer  than  350 
meters.  This  choice  appears  to  correspond  to  worst-case  guessing.  The  results  for  error 
calculation  are  obtained  for  0  degrees,  30  degrees,  and  45  degrees  at  a  detector 
height  of  5  meters.  The  actual  temperature  and  convergence  parameter  6  are  chosen  to 
be  353.54  K  and  0.001,  respectively.  The  error  of  range  determination  is  shown  in  Figure 
24,  and  the  error  of  temperature  determination  is  shown  in  Figure  25.  Both  errors  in- 
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crease  as  the  actual  range  increases.  However,  the  maximum  error  in  temperature  and 
range  is  insignificantly  small.  The  maximum  error  in  ranging  is  0.09  %  and  the  maximum 
error  in  temperature  is  0.0008  %.  Tliis  results  shows  that  the  method  performs  remark¬ 
ably  well  for  the  conditions  tested. 

In  the  calculation  of  the  transmission  factor  for  the  /-th  band,  the  precise  should 
be  calculated  according  to  the  weighted  average 


T)dk 


(5.27) 


rather  than  the  standard  average  (5.12).  It  should  be  noted  that,  the  flatter  the  photon 
emitiance  distribution  within  the  spectral  bands,  the  smaller  the  discrepancy  between 
calculations  (5.12)  and  (5.27). 

Despite  the  simplifying  assumptions  used  in  the  model,  the  results  demonstrate  the 


feasibility  of  extracting  both  target  range  and  temperature  from  thermal  radiation 
measured  over  three  bands. 


Figure  15.  Deflnition  of  vertical  line  from  the  earth  through  the  receiver 
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vertical 


Figure  17.  Deflnition  of  layers  in  atmosphere 


52 


53 


Figure  19.  Geometry  of  ranging 


Figure  20.  Geometry  of  ranging  (cosine  law) 
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Figure  22.  Graphical  solution  of  range  and  temperature 
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Figure  25.  Error  of  tempcnture  dcterminatioD 
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VI.  CONCLUSIONS 


Previous  work  on  dual  spectral  band  measurements  of  radiant  emittance  demon¬ 
strated  that,  under  somewhat  ideal  conditions,  the  ratio  of  the  measurements  could  be 
used  to  determine  the  temperature.  The  main  stipulated  condition  was  that  the  bands 
were  not  self-contained.  In  the  Chapter  III  of  this  thesis,  the  same  conclusions  were 
obtained  for  measurements  of  photon  emittance. 

A  general  model  encompassing  arbitrary  wavelength-dependent  emissivity, 
responsivity,  and  transmission  factor  was  proposed  in  Chapter  IV.  Although  the  general 
problem  is  not  analytically  tractable,  the  methods  applicable  to  constant  emissivity, 
responsivity,  and  transmission  factor  (Chapter  III)  are  approximately  valid  if  the  spec¬ 
tral  bands  are  narrow  enough.  However,  from  the  practical  view  point  of  making  a 
signal  measurement  in  the  presence  of  noise,  the  spectral  bands  should  not  be  too  nar¬ 
row.  Therefore,  the  more  general  model  is  often  needed  under  practical  conditions.  For 
reason  explained  in  Chapter  IV,  the  photon  detector  responsivity  introduced  is  fairly 
insensitive  to  wavelength  up  to  the  cut-off  wavelength.  In  addition,  for  a  broad  class 
of  materials  emissivity  curves  are  quite  fiat.  Therefore,  the  focus  of  the  reminder  of  the 
thesis  concentrated  on  effects  due  to  atmospheric  attenuation.  Specifically,  knowledge 
of  the  range  to  the  target  would  be  necessary  to  apply  the  method  discussed  in  Chapter 
IV. 

The  concept  presented  and  tested  in  Chapter  V  is  created  from  the  requirement  that 
the  dual  spectral  band  method  of  temperature  determination  would  require  range  infor¬ 
mation  to  determine  the  temperature.  This  assumes  a  disparity  in  attenuation  effects 
between  bands.  By  employing  three  bands  instead  of  two,  this  problem  is  circumvented. 
Additional  complications  created  by  changing  atmospheric  conditions  with  trajectory 
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path  are  handled  within  the  model.  The  numerical  atmospheric  attenuation  data  used  in 
the  model  testing  was  taken  from  LOWTRAN.  The  results  for  this  preliminary  model 
are  extremely  encouraging.  Future  work  needs  to  be  performed  to  develop  and  test  a 
more  rigorous  model. 
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APPENDIX  A.  DERIVATIVE  OF  ^  WITH  RESPECT  TO  T 

The  purpose  of  this  appendix  is  to  prove  equation  (3.7)  which  appears  in  the  main 
text.  Equation  (3.7)  had  been  introduced  in  {Ref.  8:  pp.  1257]  without  proof.  Since  the 
ratio  of  photon  emittance  is  a  function  of  a,  it  requires  a  few  steps  to  get  its  derivative 
with  respect  T.  For  the  reader's  convenience,  relations  (3.2a),  (3.5)  and  (3.6)  are  re¬ 
produced,  which  show  that 


fflocj 

J  gi^)d^ 


«3 


(^1) 


where 


gW 


cxp(iA)  -  1 


iA2) 


and 


a 


(^3) 


Following  the  product  rule,  the  derivative  of  the  ratio  of  photon  emittance  (Al)  with 
respect  to  temperature  is 
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'*8  flUfJl  CUfj 


where  the  temperature  dependence  is  expressed  through  a.  One  approach  to  simplUVing 
the  expression  is  to  transfer  the  dependence  of  the  integral  limit  on  a  evident  in  (A4)  to 
the  integrand.  A  substitution, 


^  ox 


iASa) 


will  transfer  this  dependence.  It  follows  that 


f  ^g{ax)dx. 


Therefore  according  to  the  product  rule, 


■  J  ^gWd^l^  -  J  ^£{ax)dx  +  «  |  g(ax)rfx. 

•'**(  •’jc,  J 


Now  from  (A2)  and  the  transformation  rule  (A5a), 


cxp(ax)  -  1  ’ 


and  it  follows  from  (A7)  that 
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Jl  (  \  exp(ajf) 

exp(«x)-  1  "  [  exp(oLc)-  l]*  ’ 

Direct  8uV)stitution  of  (A8)  into  (A6)  leads  to  : 


m 


I  ^g{cuc)dx  + 


2(»)* 
exp(ajf)  —  1 


dx  — 


(ttx)^exp(ttjr) 

J  Cexp(ajr)-- Ip 


dx.  (/49) 


Noting  the  equivalence  between  g(ouf)  defined  in  (A7),  and  the  integrand  of  the  second 
integral  (A9),  it  follows  that : 


exp(<xjr) 

C  cxp(ax)  - 1]^ 


dx. 


(AlO) 


After  reversing  the  transformation  (A5a),  fhe  integiands  in  (AlO)  will  depend  on  ij/  and 
(AlO)  leads  to: 


•’oW/ 


/*«/ 


’'ax, 


exp(^) 

C  exp((/')-  1]^ 


dilf. 


iAU) 


For  compactness  of  representation,  define 


A^) 


fj/^  exp(</>) 
C  exp(^/^)  - 


(A  1 2) 
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Then  combining  the  result  (A3)  with  (A1 1)  leads  to  : 


.  ■  -L I  (/<13) 

L  J 


After  substituting  (A  13)  into  (A4) 


g{^)d4> 


(-^14) 


which  after  simplification  leads  directly  to  (A16)  or  (3.7),  as  given  in  the  text. 


dT  T 


PeLXj  rax,  faxj  Pax, 

M)d\J/  g{4f)dflf  “  A^)d^ 

1  ‘[a _ _ 

7'  pax, 

gW^ 

.  “I 


(^15) 
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APPENDIX  B.  ASYMPTOTIC  APPROXIMATION  FOR  RATIO 

CALCULATION 

The  purpose  of  this  appendix  is  to  prove  equations  (3.22)  and  (3.23)  which  appear 
in  the  main  text  where  the  asymptotic  approximations  are  applied. 

A.  FOR  LARGE  ALPHA 

For  large  a  (i.e.,  small  temperature),  since  x,  values  are  fixed,  0  takes  on  large  posi¬ 
tive  values  Therefore,  functions  are  approximated  from  (3.5)  and  (3.8)  as 

follows: 


cxp((A) 


(51) 


and 


cxp(!/')  ’ 


After  direct  substitution  into  (3.8),  it  follows  that: 


(52) 


f,g]^ 


1  cxp(  -  cxp(  - 


naxj  /•(«« 

exp(  -  ^)#  exp(  - 


(53) 


One  round  of  integration  by  parts  yields 


-j*  ^A*exp(-^^y^j[-^’exp(-^)]^  +  3j  exp(  -  j 


After  algebraic  simplification,  (B4)  becomes 


paJC, 

exp(  -  exp(  ~ 

*'«x, 

f«a 

cxp(  ~  ~  exp(  - 


(55) 


According  to  (Ref.  12:  pp.  85] 


f 


exp(  -  \lf)dili  -  -  exp(  -  ^)  -  exp(  -\lf)-2  exp(  -  0),  (56) 


which,  for  large  positive  0  (B6),  is  approximately  given  by 


0^  exp(  -  0)(/0»  -  0^  exp(  -  0). 


(57) 


After  substitution  of  the  approximate  form  (B7)  into  the  integrals  (B5),  [f,g]  becomes 
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exp(  -  exp(  -  exp(  -  exp(  - 


which,  aAer  simplification,  becomes 

l/',gl“C(aJf2)^  exp(  -  a^fj)  -  exp(  -  a)3C(ajf4)^  cxp(  -  01X4)  -  exp(  -  ouj)] 

-  [(axj)^  exp(  -  axj)  -  exp(  -  a)3C(«X4)^  exp(  -  a;f4)  -  exp(  -  axj)] 

Since,  for  large  a,  the  second  terms  within  each  bracket  of  (B9)  become  dominant,  [f,g] 
is  approximated  as 

(/;g]fea®[x3  exp(  -  0LX3)  exp(  -  a)  -jf|  exp(  -  ajfj)  exp(  -  «)].  (510) 

Recombination  of  (Bll)  leads  to: 

-  X3)  exp  I  -  a(jf3  +1)].  (511) 

This  last  form  is  applied  in  Section  C  of  Chapter  III. 

B.  FOR  SMALL  ALPHA 

For  small  a  (i.e.,  large  temperature  ),  \j/  is  small.  Thus,  gi\j/)  and  are  approxi¬ 
mated  from  (3.5)  and  (3.8)  by  binomial  expansion.  After  application  of  binomial  expan¬ 
sion,  g{^)  is  integrated  as 


(BU) 


And  a  binomial  approximation  for  small  ^  leads  to: 

jgWd^^j  1^(1  -  -Y 


(513) 
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which  is  integrated  as 


(B14) 


Similarly,  after  applying  binomial  expansion  to  (3.8),  the  integration  becomes 


(515) 


Algebraic  expansion  of  (B15)  leads  to: 


(516) 


which  is  approximated  as 


2  16 


(517) 


Therefore,  for  small  oc,  [f,g)  becomes 


(518) 


which  is  expanded  as 
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UtSJ- 


—  ^3/ —  22  V'J  “*3/ 


V'2  *  ^v-4  "i  / 


-  -j^  (xi  -  xixx!  -1) + -^  «‘(4  -4x^1  -')-■«(***-  ^?x4  - ') 


(819) 


Since,  for  small  alpha,  a’  terms  are  dominant,  (f,g]  approximated  as 

l/;^l=^  4"  -i)  -  (^?  -  l)(^i  -  4)]-  .  (520) 


This  final  form  is  applied  in  Section  C  of  Chapter  III. 


APPENDIX  D.  ALGORITHM  TO  CALCULATE  THE  RATIO  OF 
OUTPUT  VOLTAGE  AT  THE  DETECTOR 
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APPENDIX  E.  MODEL  FOR  PASSIVE  RANGE  AND  TEMPERATURE 

DETERMINATION 

A.  MAIN  ALGORITHM 


calculation 

temperature 


Ideal  ratio 

8eeEQ(5.7) 


B.  LOWTRAN  DATA  GENERATION 


Inputs 
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C.  LAYER  (SUBROUTINE) 


Inputs 
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D.  MEASURED  RATIO  PREDICTOR 


Inputs _ 

Range  R1  ,R2 
Temperature  T 
Angle 
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